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We analyze superconductivity in the cuprates near the onset of an incommensurate charge den¬ 
sity wave (CDW) order with momentum Q = (Q, 0)/(0, Q), as observed in the experiments. We 
first consider a semi-phenomenological charge-fermion model in which hot fermions, separated by 
Q, attract each other by exchanging soft CDW fluctuations. We find that in a quantum-critical 
region near CDW transition, T c = Ag c , where g c is charge-fermion coupling and A is the prefactor 
which we explicitly compute. We then consider the particular microscopic scenario in which CDW 
order parameter emerges as a composite held made out of primary spin-density-wave fields. We 
show that charge-fermion coupling g c is of order of spin-fermion coupling g s . As the consequence, 
superconducting T c is substantially enhanced near the onset of CDW order. Finally we analyze the 
effect of an external magnetic held H. We show that, as H increases, optimal T c decreases and the 
superconducting dome becomes progressively more conhned to the CDW quantum-critical point. 

These results are consistent with the experiments. 


I. INTRODUCTION 

Understanding of the nature of charge order in high- 
T c cuprates and of its effect on superconductivity is es¬ 
sential for the full understanding of the complex physics 
in this materials. Charge order has been observed in 
the cuprates quite some time ago^, but was originally 
though to be present only in La-based materials. Recent 
wave of discoveries of an incommensurate charge-density- 
wave (CDW) order in Y-, Bi-, and Hg- based cuprates^— 
has demonstrated that charge order is ubiquitous to all 
families of high-T c cuprates. A true long-range CDW or¬ 
der has so far been observed only in a finite magnetic 
field, but short-range static order (probably pinned by 
impurities) exists already in zero field. On a phase dia¬ 
gram, CDW order been detected within the pseudogap 
region, and its onset temperature is the highest around 
doping level x ~ 0.12. The CDW has an incommensu¬ 
rate momentum Q = Q y = (0, Q) or Q x = ( Q , 0) and the 
order is likely uni-axial, i.e., it develops, within a given 
domain, with only Q x or Q y (Ref. 0 ). 

The discovery of the CDW order raised a number of 
questions about its originiir— , a potential discrete sym¬ 
metry breaking before a true CDW order sets mi i 14 : 22 " — , 
and the relation between CDW order (or its fluctuations) 
to pseudogap behavio r 14 ! 18 : 20 ! 27 and Fermi surface (FS) 
reconstruction— 

In this paper we discuss another issue related to CDW 
its effect on superconductivity. We take as inputs three 
experimental facts. First, T c , as a function of doping, 
has a dip or a plateau at around x ~ 0.12, where the 
onset temperature of CDW is the larges t 30 ! 31 . Second, 
when CDW is suppressed by applying pressure^, super¬ 
conducting T c increases. Third, when a magnetic field is 
applied^, the dip grows and at high enough field the su¬ 
perconducting dome splits into two, and the one at larger 


dopings is centered at the same x at which CDW order 
develops at T = 0 (Ref. [H). In other words, supercon¬ 
ductivity forms a dome on top of quantum-critical point 
(QCP) for the onset of the CDW order. 

The first two results can be naturally understood if we 
assume that CDW and d -wave superconductivity are just 
competing orders, i.e., when one order is at its peak, the 
other one is suppressed. The third observation, on the 
other hand, requires one to go beyond a simple “com¬ 
peting order” scenario because the presence of the dome 
of superconductivity on top of CDW QCP indicates that 
superconductivity is at least partly caused by soft fluctu¬ 
ations of CDW order which then must develop at larger 
energies than the ones related to superconductivity. 

In our analysis we explore CDW-mediated supercon¬ 
ductivity in some detail. We perform our analysis in 
two stages. At the first stage we put aside the issue 
what causes CDW order, assume that this order devel¬ 
ops below some critical doping, and consider a semi- 
phenomenological model of fermions interacting by ex¬ 
changing soft CDW fluctuations with momenta Q. This 
model is quite similar to the spin-fermion model, consid¬ 
ered in earlier studies of spin-mediated superconductiv¬ 
ity for cuprates, Fe-pnictides, and other correlated mate¬ 
rials^—, and we dub this model the “charge-fermion 
model”. The charge-fermion and spin-fermion models 
are similar but differ in detail because of the difference 
between the CDW momentum Q and antiferromagnetic 
momentum (7r,7r), and also because of the difference in 
the spin structures of charge-mediated and spin-mediated 
interactions (spin Kronecker 5 functions vs spin Pauli ma¬ 
trices). One qualitative consequence of these differences 
is that charge-mediated interaction gives rise to super¬ 
conductivity with a full gap in each region where the FS 
crosses the Brillouin zone boundary (anti-nodal regions 
in the cuprate terminology), but it does not couple su¬ 
perconducting order parameters from different regions, 
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hence it alone cannot distinguish between s-wave and d- 
wave pairing symmetries— 

We perform quantitative analysis of the pairing 
within the charge-fermion model in the most interesting 
quantum-critical regime right above QCP for CDW or¬ 
der. In this regime, charge-mediated interaction gives 
rise to the pairing but also destroys coherence of hot 
fermions (the ones at the FS separated by Q). Super¬ 
conducting T£ h is determined by the interplay between 
strong tendency towards pairing and strong pair-breaking 
effects associated with the self-energy. We compute Lan¬ 
dau damping of soft bosons and fermionic self-energy, 
and then obtain and solve the linearized gap equation 
with renormalized fermionic and bosonic propagators. 
We show that in the quantum-critical region T c is finite 
and scales with the effective charge-fermion coupling con¬ 
stant g c : T/ h = A c g c , where A c « 0.0025. 

At the second stage we consider the specific micro¬ 
scopic scenario for CDW order - the one in which CDW 
is a composite order parameter made out of primary 
( 7 r, 7 r) spin fluctuations. Within this scenario, spin fluc¬ 
tuations are assumed to develop first, at energies com¬ 
parable to bandwidth, while CDW fluctuations develop 
at smaller energies and do not provide substantial feed¬ 
back on spin fluctuations. This composite order sce¬ 
nario requires fermion-fermion interaction to be compa¬ 
rable to the bandwidth (otherwise spin fluctuations do 
not develop at high energies) and inevitably is partly 
phenomenological. We will not discuss a complemen¬ 
tary, renormalization group-based, truly weak coupling 
scenario in which all fluctuations (spin, charge, super¬ 
conducting) develop simultaneously at low energy and 
mutually affect each other—. 

Various versions of magnetically induced charge bond 
and site orders have been proposed over the last few 
yearsiirii^-is^i^^ some focused on CDW with di¬ 
agonal momentum ( Q, Q ), and others on CDW with mo¬ 
menta Q x = (Q,0) and Q y = (0, Q). Motivated by ex¬ 
periments, we base our discussion on the soft fluctuations 
of axial CDW with momenta Q x or Q y . Previous studies 
have founds that axial CDW has a partner - an incom¬ 
mensurate pair-density wave (PDW), and fluctuations in 
CDW and PDW channels develop simultaneously. To 
keep presentation focused, we concentrate on CDW and 
neglect PDW fluctuations. The latter can in principle 
also mediate pairing interaction, but are unlikely to de¬ 
structively interfere with CDW fluctuations. 

In an earlier work— we have shown that the axial CDW 
order develops in a paramagnet at a finite T c d w , provided 
that magnetic correlation length £ s is large enough. As 
£ s decreases, T c( j w also decreases and vanishes at some 
finite £ s ,cr) setting up a CDW QCP at a finite distance 
from a magnetic QCP (the one at = 00 ). Near 6,,cr 
CDW fluctuations become soft and give rise to singular 
pairing interaction mediated by these fluctuations. We 
show the behavior of T c d w vs £ s schematically in Fig. |T] 

Whether this additional pairing interaction substan¬ 
tially affects superconducting T c depends on the inter- 



FIG. 1. Doping range where CDW order with momentum 
( Q , 0) or (0, Q) emerges within itinerant spin-fluctuation sce¬ 
nario. The effective interaction in the CDW channel is made 
out of two spin-fluctuation propagators. The CDW order de¬ 
velops only when the system is sufficiently close to a magnetic 
instability and terminates at doping x? r ‘. different from a;®? for 
antiferromagnetism. In the shaded region near localiza¬ 
tion of electronic states (Mott physics) becomes relevant (Ref. 
[45l l. and spin-fluctuation approach needs to be modified. In 
this region, the onset temperature of CDW order drops. 


play between charge-fermion coupling g c and the cou¬ 
pling g s between fermions and primary spin fluctuations. 
The argument is that spin-fluctuation exchange by it¬ 
self gives rise to superconducting pairing, and at large 
£ s the corresponding T/ p scales with g s : T/ p = A s g Sl 
where A s « 0.007 (Ref. H^). As £ s decreases, T® p also 
decreases but remains finite. The two-dome structure of 
T c (x), observed in the cuprates in an applied magnetic 
field, can be understood within this approach only if near 
= £s,cr, Tc h — T°p. If this does not hold, i.e., Tj? h is 
smaller than T/ p , the contribution to superconductivity 
from charge-mediated exchange is only subdominant to 
the one from spin fluctuations. In this situation, the 
only effect on T c from CDW is due to direct competition 
between CDW and superconducting orders. This com¬ 
petition may give rise to additional reduction of T c in 
a magnetic field, given the experimental evidence that 
CDW order increases in a field. However, it cannot give 
rise to a peak of T c above CDW QCP. 

We evaluate charge-fermion coupling g c within the 
RPA-type analysis of charge fluctuations near CDW 
QCP: uf{q) cx £4/(1 — U c Il c (q)) = g c / (£“ 2 + (9 — Q) 2 ), 
and show that g c is comparable to spin-fermion coupling 
g s . This result may look strange because charge fluc¬ 
tuations, viewed as composite objects made out of spin 
fluctuations, develop only in a narrow range near the FS 
points separated by Q x or Q yi with the width in mo¬ 
mentum space of order A ~ ^T 1 . As the consequence, 
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g c ~ t4C7 2 - However, the “bare” interaction in the 
charge channel, U c , is a composite object made out of two 
spin fluctuation propagators and two fermionic Green’s 
functions (see Fig. |SJ). This composite object behaves as 
g s £ 2 . As the consequence, U c ^~ 2 is not reduced by 
and g c differs from g s only by a numerical factor. 

To properly calculate the ratio g c /gs one needs to do 
full-scale dynamical calculations, even if one restricts 
with ladder series of diagrams, like in RPA. In this work, 
we use a simplification and approximate the bare inter¬ 
action in the charge channel U c by a constant within 
the momentum range A ~ around proper Fermi sur¬ 
face points (hot spots) and set it to zero outside this 
range. We compute the polarization operator n c (g, fl m ) 
and use RPA to obtain charge-mediated effective inter¬ 
action within fermions. We use the condition for CDW 
QCP: 1 = U C H C (Q) to fix A and obtain explicit rela¬ 
tion between g c and g s . Within this approach, we find 
ft ^ Js- 

We argue, based on this analysis, that the enhance¬ 
ment of superconductivity near a CDW QCP is sub¬ 
stantial, i.e superconductivity in the cuprates comes 
from both spin and charge fluctuations. We present the 
schematic phase diagram in Fig. [2j This scenario explains 
the developments of two domes of T c in a high magnetic 
field: one, at smaller doping, is due to critical spin fluc¬ 
tuations, and another, at larger doping, is due to critical 
charge fluctuations. 


1. Relation to earlier works 

The pairing by charge fluctuations has been studied 
before. In the context of the cuprates, DiCastro, Castel- 
lani, Grilli, and their collaborators^ analyzed in detail 
the pairing mediated by axial CDW fluctuations near the 
onset of charge order, which was assumed to develop on 
its own rather than being induced by SDW fluctuations 
They found that charge-mediated 4-fermion interaction 
is attractive in both d -wave and s-wave channel and does 
not distinguish the two. They argued that some other 
mechanism, e.g., antiferromagnetic spin fluctuations, lifts 
the degeneracy in favor of d-wave. We obtain the same 
results in Sec. II. The novel part of our analysis in this 
Section is the calculation of charge-mediated T c in the 
quantum-critical regime. 

This pairing problem near CDW QCP has certain sim¬ 
ilarities to the pairing at the onset of a nematic order, 
which has been extensively studied in recent years^Zr— . 
It has been argued that Q = 0 nematic fluctuations en¬ 
hance all partial components of the pairing susceptibility. 
The case of QCP at small Q x /Q y is less unrestrictive in 
this respect but still, s-wave and d -wave channels are de¬ 
generate for CDW fluctuations. 

The analysis of the electron-mediated pairing near a 
QCP for density-wave order is also quite interesting from 
a general perspective as it adds one important new el¬ 
ement not present for the pairing away from a QCP. 



FIG. 2. Schematic phase diagram for the interplay between 
superconductivity mediated by spin and charge fluctuations. 
Panel (a): the onset temperature T’jP of spin-mediated super¬ 
conductivity, as a function of doping. We assume that im¬ 
purities kill superconductivity above a certain doping. Panel 
(b): the onset temperature T^ p of charge-mediated supercon¬ 
ductivity near the onset of CDW order. We compute 
in this work. Panel (c): The full phase diagram. The non¬ 
monotonic T c (x) is obtained by combining spin-mediated and 
charge-mediated contributions to T c from panels (a) and (b) 
(dashed lines). 


Namely, the same interaction which gives rise to strong 
attraction also destroys fermionic coherence and prevent 
fermions from developing supercurren t 51 ! 52 . Supercon¬ 
ductivity then may or may not emerge depending on the 
interplay between these two opposite tendencie s 47 ! 53 . 

The rest of this paper is organized as follows. In Sec. II 
we introduce and analyze charge-fermion model of itiner¬ 
ant electrons coupled to near-critical CDW fluctuations. 
In Sec. Ha we derive bosonic and fermionic self-energies. 
In Sec. lib we study the pairing problem and obtain T“ h 
due to charge-fluctuation exchange near a CDW QCP. 
We show that Tf 1 scales with the charge-fermion coupling 
constant g c . In Sec. Ill, we relate g c and spin-fermion 
coupling g s within the magnetic scenario for CDW. In 
this scenario, a CDW order parameter field emerges as a 
composite object made out of two spin-fluctuation prop¬ 
agators. We show that g c is comparable to g s and may 
even exceed it. In Sec. IV we discuss in some more de¬ 
tail superconducting dome around CDW QCP. Finally in 
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FIG. 3. The Fermi surface of a 2D electronic system on a 
square lattice and location of hot spots for charge-mediated 
interaction. CDW “hot” spots are defined as points on the 
Fermi surface separated by either Q y = (0, Q) or by Q x = 
( Q , 0). In a generic situation, there are 8 hot spots for Q y and 
eight for Q x . Motivated by experiments, we consider the case 
when hot spots for Q y and for Q x merge. In this situation, 
CDW hot spots coincide with hot spots for (7r, 7r) magnetism, 
and there are eight of them on the Fermi surface We label 
these hot spots by 1-8 (5 = —2, 6 = —1, etc). The arrow 
shows the direction of the Fermi velocity at hot spot 1. Fermi 
velocities at other seven hot spots are related to this one by 
symmetry. Near CDW instability hot fermions interact with 
each other by exchanging CDW fluctuations with momenta 
Qx, y - Here and in other figures below we represent charge 
fluctuations by dashed lines. 


Sec. V we discuss the results and present our conclusions. 


II. THE CHARGE-FERMION MODEL 

We begin with a semi-phenomenological analysis. We 
assume, without specifying the reason, that CDW order 
with momentum Q y = (0, Q) and/or Q x = (Q, 0), devel¬ 
ops at some critical doping x™, and consider the model 
of 2D fermions interacting by exchanging near-critical, 
soft charge fluctuations. We dub this model as “charge- 
fermion model” by analogy with the spin-fermion model 
which was introduced to describe a system near a mag¬ 
netic QCP. As our goal is to describe low-energy physics 
(energies well below the bandwidth), we focus on mo¬ 
mentum ranges around the CDW “hot” spots on the FS, 
defined as points which are separated by CDW momen¬ 
tum Q x or Q y . For a generic Q there are 16 CDW hot 
spots, eight corresponding to Q y and another eight cor¬ 


responding to Q x . For simplicity we consider the case 
when Q is such that CDW hot spots are at or near the 
crossing between the FS and symmetry lines in the Bril- 
louin zone k x ±k y = ±7r. Then hot spots from Q y sector 
merge with hot spots from Q x sector, and the total num¬ 
ber of hot spots become eight. We label these points as 
1-8 in Fig. [3] This approximation works reasonably well 
for the values of Q extracted from the experiments^. 

We linearize the fermionic dispersion in the vicinity 
of a hot spot i as e i j. = v* • kj, where v, is the Fermi 
velocity and is the momentum deviation from the hot 
spot i. We define the Fermi velocity at the hot spot 1 as 
v i = (%,%), the velocities at other hot spots are related 
by symmetry. In the cuprates, the FS at the antinodal 
region is “flattened”, and we have v x < v y (see Fig. [3j. 
We define a = v x /v y < 1. 

The action of the charge-fermion model can be written 
as 

S = J dk^24 a (k)(-iu m + e iX .)ci ta (k) 

i,a 

+ I dqxociv) Y 

a—x,y 

+ 9c Y fy [ dkdk'4oc(k) c i+i, a (k')<pl(tf - k) 

i= 1,3,5,7; ^ 

ol 

+ 9c Y fx f dkdk' 4 a &)ci+^a{k')(4 x {k - k') 

2 = 1 , 2 , 3 , 4 ; ** 

a 

+ h.C; (1) 

where Ci a is a fermion field with i labeling hot spots 
and a labeling spin projections. Hot spots i and i + 1 
are separated by CDW momentum Q yi and hot spots 
i and i + 4 are separated by Q x (see Fig. O. The 
scalar field <p XtV is a CDW order parameter field with 
momenta near Q x /Q y . In Eq. dTJ) we have used short¬ 
hands k = (ui m , k) and q = (f2 m ,q), where cj m (fl m ) are 
fermionic (bosonic) Matsubara frequencies. The bosonic 
momentum q is measured as the deviation from CDW 
momenta Q x or Q y , and the fermionic momentum k is 
measured as the deviation from the corresponding hot 
spot. The form-factors f x y determine relative magnitude 
and sign between CDW orders in different hot regions. 
In general, the CDW order has both d-wave and s-wave 
components. A pure d -wave order would correspond to 
fy ’ 5 = - fy ’ 7 = 1 and /x 1 ’ 2 = -/* 3 ’ 4 = 1 - 

We assume, like it is done in the spin-fermion model, 
that static charge susceptibility comes from fermions 
with energies larger than the one relevant to super¬ 
conductivity and approximate it by a simple Ornstein- 
Zernike form Xc = X 0 c/{ql + q% + £)T 2 ); where £ c is the 
CDW correlation length. The prefactors for q x and q y 
may in general differ— because Q x and Q y are not along 
Brillouin zone diagonals, but this difference can be ab¬ 
sorbed into rescaling of q and does not affect our analysis. 
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The coupling g c is a phenomenological charge-fermion 
coupling constant. The corresponding charge-mediated 
4-fermion interaction term in the Hamiltonian is (we set 
Q = Q y for definiteness) 


Td'vat — U c (q) 'y ' Cfc ^Cp ^Ck—q^sCk+q’pSapfi'yS (2) 


with 


Uf{q) =9cXc{q) = — 


9c 


£c + (q — Qy ) 2 


( 3 ) 


The effective coupling g c = g^Xoc, and the sign conven¬ 
tion is such that the interaction appears with a factor — 1 
in the diagrammatic theory. 

The charge-fermion model is defined self-consistently 
when its fluctuations cannot modify the physics at lat¬ 
tice energies, and this requires that g c must be small 
compared to the fermionic bandwidth. 


A. Normal state analysis 


i 



i + 1 


(a) 


Qy 



l 



l+l 




i + 4 


i 


(c) 


FIG. 4. The one-loop diagrams for the bosonic self-energy 
tlx,y [Panels (a), (b)] and fermionic self-energy E i [Panel (c)[. 


For the computation of superconducting T)r h , mediated 
by charge fluctuations, we need to know normal state 
properties. We use Eq. m and compute self-energies 
for the bosonic fields (f> x and <p v and for the fermionic 
field We show corresponding diagrams in Fig. [I] We 
compute basonic and fermionic self-energies in a self- 
consistent fashion, like it was done in the earlier works 
on the spin-fermion mode l 11 ’ 52 ' 55 . Namely, we first eval¬ 
uate one-loop bosonic self-energy (the bosonic polariza¬ 
tion operator) using free fermions and show that it has 


the form of Landau damping, then use the full dynam¬ 
ical bosonic propagator to calculate one-loop fermionic 
self-energy and show that it is strong but predominantly 
depends on frequency, and then verify that frequency de¬ 
pendent fermionic self-energy does not affect the Landau 
damping. This self-consistent procedure becomes exact 
if we neglect subleading terms in the self-energy, which 
depend on the fermionic dispersion e&. This can be rig¬ 
orously justified if we extend the model to M fermionic 
flavors and take the limit M —I oo (e.g. Ref. EH) , or ex¬ 
tend the number of pairs of hot spots from 4 to IV and 
take the limit IV —> oo (e.g. Ref. fS2h . We will use the 
latter extension to justify our analysis. 


1. bosonic polarization operators 


The bosonic polarization operator for cj) y is given by 
the diagram in Fig. [Ha), and the expression for cf> x is 
related by symmetry. The full polarization operator is 
a sum of contributions from four pairs of hot fermions 
which are separated by Q y . These pairs are (1,2), (3,4), 
(5,6), and (7,8). There are no Umklapp processes for 
incommensurate CDW order, in distinction to SDW case 
in which Q = ( 7 r, 7 r) and Umklapp processes are allowed. 

For the contribution to n from fermion pair (1,2) we 
obtain 


ni(n m , q) — —2<? 2 T Gi(w m , k)G 2 (w m + fl m , k + q) 

( 4 ) 


where the factor 2 comes from summation over spin in¬ 
dices. To simplify the notations we will drop the tildes 
on momenta hereafter. The Green’s functions are given 
by Gi 5 2 = 1 — £ 1 , 2 ) and fermionic dispersions can 

be written as ei^ = vi^ • k. We transform the mo¬ 
mentum integral dk x dk y into de\dc 2 by adding the Ja¬ 
cobian J = l/|vi x V 2 I = 1 /vp x (a 2 + l)/(2a), where 

y F = ^J v x + . Because this Jacobian is independent 

of q, the polarization operator in this approximation is a 
function of frequency only. We subtract from n 1 (U m ) its 
frequency-independent part ni(0), which only renormal¬ 
izes the position of the CDW QCP and is not of interest 
to us. The subtraction makes integral over internal fre¬ 
quency convergent, and evaluating the integrals we ob¬ 
tain at zero temperature 
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IIi(f2 m ) — — 


9c ot 2 + 1 


47T 3 U^ 

a J 

9c « 

2 + l f 

dnvp 

a J 

9c 

a 2 + ! 


dLo m deide2 


Airv^, a 


(iu> m - ei )[i{u> m + £l m ) ~ £ 2 ] (*w m - ei )(iu) m - e 2 ) 

diO m [sgn (u m ) sgn (u m + n m ) - 1 ] 


(5) 


This is a conventional Landau damping term. We do 
the same calculation for hot spot pairs (3,4), (5,6), and 
(7,8). Because the Jacobians l/|vi x v 2 | = l/|v 3 x V 4 I = 
l/|v 5 x v 6 1 = l/|v 7 x Vg| = 1 /v 2 F x (a 2 + l)/( 2 a), all 
contributions are the same as ([5]). Therefore n y (fl m ) = 
4IIi(n m ). It is easy to verify that for <fi x the self-energy 
is the same as for (jj y , i.e., H x = Tl y . 

Including the polarization operator into the propaga¬ 
tors for (f) x and (j) y fields, we obtain 


Xc(^rm q) 


_ XOc _ 

£c 2 + 9 X + gy + 7e|^m| 


( 6 ) 


where 


4 g c a 2 + 1 
47 Wp a 


(7) 


and we recall that g c = g 2 X 0 c- The overall factor of 4 
in the numerator of 0 is the number of pairs of hot 
fermions. To extend the model to large N one just has 
to replace 4 by N. We will use this extension below. 

The total dynamical interaction mediated by charge 
fluctuations can then be written as 


I7 c eff (q,^m) =g 2 c x( q) 


_ 9c _ 

£c 2 + (q — Qz) 2 + 7c|flm| 

( 8 ) 


2. fermionic self-energy 


We now use Eq. 0 and evaluate one-loop fermionic 
self-energy. The corresponding diagram is presented in 
Fig. 0c). Observe that for any hot fermion, interactions 
mediated by both bosonic fields, <f> x and </> y , contribute 
to the self-energy. For example, for a fermion at hot spot 
1, 4> x and 4> y scatter it to hot spots 2 and 5, respectively. 

The self-energy depends on the location of a fermion on 
the FS and on the distance to CDW QCP. Below we will 
be interested in superconductivity right at CDW QCP, 
hence we will need the self-energy right at this point. 
Accordingly, we set = 0 in the charge fluctuation 
propagator. 

For the self-energy contribution to a fermion at hot 


spot 1 from Q y scattering, we have 


E y (k ,u m ) = T J2 C^f(k-k , ,w TO -OG 2 «,k , ) J 

(9) 


where k is the deviation from hot spot 1. We place k on 
the FS, i.e., set k± = vy • k = 0, which gives k y = — ak x . 
At T = 0, we rewrite Eq. 0 as 


Ey 0 ||, co m ) 


-gc 

87T 3 


0|| - fcy) 2 + O-L 


duj^dk^dk^ | 

iuj' m — vpk' ± 

1 

- k' ± ) 2 +7c |UJm 



( 10 ) 


where k' ± and fc|, are perpendicular and parallel compo¬ 
nents of k' with respect to the Fermi surface at hot spot 
2, i.e., k' ± = v 2 • k', and fc|| and k± are components of 
external k, defined relative to the FS at the hot spot 2, 
i.e., kj_ = v 2 • k, where v 2 = (v x , —v y ). Using k y = —ak x 
we obtain 


k± = 2ak\\/(a 2 + 1). (11) 

We integrate over k' ± first, and complete the integra¬ 
tion contour over the half plane with only one pole. We 
obtain 


E. y (fc||, UJ m ) 



0|| - |) 2 +7c|w m 


dk\\du' m sgn(w( ra ) 

\J{h - k' n ) 2 +7 c |w m - uj' m \ 

1 

U'm\ + [\u'm\l v F + tk± Sgn(w^)] 

( 12 ) 


We will see that typical internal frequencies u>' m are of 
the same order as external uj m . Then, at small enough 
frequencies one can safely neglect the \ui' m \/vF term in 
the denominator. With this simplification we obtain 


Ey (fc||, UJm) 

= ig c f _ 

8tt 2 vf J 0|| - / 

= 9 sgn(w m ) 

2'KVF'fc 


dk'\\ du}' m sgn(w( n ) 

:') 2 + 7 c | w m + fc 2 

\jxc\^m\ +k\— | fc_L I 


(13) 
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Plugging Eq. m into Eq. m we finally obtain 


Ej / (fci|,w m ) = %9c —sgn(w m ) 

ZnvF'y c 


Hc\Urr 

2 wf a 
TV a 2 + 1 


2afc|| 
a 2 + 1 


2ak\\ 


a 2 + 1 


sgn(w m ) 




2afc|| 
a 2 + 1 


2afcii 


a 2 + 1 


(14) 


The self-energy from Q x scattering is obtained in the 
same way: 


S x (k,CJ m ) g c T 'y ^ Xc{^m ^ )Gr, 

(15) 


As Fermi velocities at hot spot 5 and 2 are antiparallel, 
we have Gs^^k') = G 2 (w(„,— k'). Comparing Eqs. 
m and (usd, we then immediately find that E^ = E y . 
Combining the two we obtain 


E(fcii,w m ) = %9c sgn(u > m ) 
fvf^c 



2ak\\ 
a 2 + 1 


4wf ol 
a 2 + 1 


sgn(a; m ) 


x 



2afcy 
a 2 + 1 

(16) 


It is easy to verify that this results holds for around all 
hot regions 1-8, and in each region fcii is the deviation 
from the corresponding hot spot along the FS. The func¬ 
tional form of the self-energy as in Eq. (1161) was first ob¬ 
tained for the spin-fermion model in Ref. |52| (for a = 1) 
and Ref. [ll] (for arbitrary a). Right at a hot spot, the 
fermionic self-energy has a non-Fermi liquid (NFL) form: 

E(0,w m ) = *sgn(w m )\/wo c |u; m |, (17) 

where ujq c = (4/ N)g c /ir x a/(a 2 + 1). Away from a hot 
spot, at fc 2 T§> 7 c|w m |, the self-energy E(w m , fc||) retains a 
Fermi liquid (FL) form at the smallest w m , i.e., we have 

s ^=w{£- r ^) + °^ <18 > 


One can easily verify that the inclusion of fermionic self¬ 
energy E(fc,w m ) oc isgn(w m ) will not change the polar¬ 
ization operator, i.e., II(f2 m ) retains the same form even 
if we compute it using dressed fermions. 

To verify self-consistency of the calculations, we also 
computed the self-energy away from the FS. We found 
E(A;_l, 0) oc ( 1/N)vpk± log(A/|np’fcj_|), where A is the up¬ 
per cutoff in momentum integration. The presence of the 
logarithm implies that Fermi velocity also acquires sin¬ 
gular renormalization at CDW QCP (Refs, [ll] and [52l) . 
This singularity breaks self-consistency of one-loop cal¬ 
culation of fermionic and bosonic self-energies if we keep 
N finite, but self-consistent procedure still remains rigor¬ 
ously justified at this loop order if we set TV —y 00 . The 
situation gets more complex at higher loop orders due 
to special role of forward scattering and backscattering 
processes which give rise to the dependence of E on k± 
without the factor 1/TV (Refe. lifil l5f)l. [Hfil. and[H3). How 
important are these effects in unclear and in this work 
we restrict with one-loop self-energy. 


B. The pairing problem 


We now use the normal state results as inputs for the 
analysis of the pairing mediated by CDW fluctuations 
with momenta around Q x and Q y . To leading order in 
1/TV, the pairing problem can be analyzed without vertex 
corrections, by summing up the ladder series of diagrams 
in the particle-particle channe l 11 : 52 . 



FIG. 5. Diagrammatic representation of the coupled ladder 
equations for superconducting order parameters $1 and $ 2 , 
which involve fermions in hot regions 1 and 2 in Fig. 1. 


We first focus on the momentum region where hot 
spots 1,2, and 5 = —2, and 6 = —1 are located (see 
Fig. [31). We introduce superconducting order parameters 
$1 ~ (cic 6 ) = (cic_i) and 4> 2 ~ (c 2 c 5 ) = (c 2 c_ 2 ) and 
obtain in the standard way a set of coupled gap equations 
for the two condensates. The interactions with momen¬ 
tum transfer Q x = ( Q , 0) and Q y = (0, Q) connect hot 
spots 1,5 with 2,6, and hot spots 1,2 with 5,6, corre¬ 
spondingly. As the consequence, the interactions relate 
<!>i with $ 2 and vise versa. We show the equation for $1 
diagrammatically in Fig. [5] In analytical form, we have 






































(19) 


$i(fc) =TY,Uf(k - k') X [G 2 {k')G 5 {-k')$ 2 {k') + G 2 {-k')G b {k')<5> 2 {-k')} 

k' 


where k = (w m ,k) and k' = (u}' m , k') and k,k' are mo¬ 
mentum deviations from the corresponding hot spots. 
The fermionic Green’s function are given by Gi(k ) = 
1 /{iu> m — £i(k) + Ej(/c)]. The equation for 4>2 in terms 
of $1 has the same form, and thus 4>i and $2 have the 
same magnitude. 

Because the two kernels in the Eq. (ITUT) (the prefactors 
for $2 hr the right hand side) are both positive (we recall 
that U e g is positive), the U(l) order parameters $1 and 
$2 have the same phase, i.e., $1 = $ 2 - By the same 
token, the SC order parameters in the momentum range 
near hot spots 3,4,7, 8 , namely 4>3 ~ {c^cs) = ( 030 - 3 ) 
and $4 ~ ( 0407 ) = ( 040 - 4 ) are also equal. The kernels 
of the gap equations in the regions 1,2,5 ,6 and 3,4,7 ,8 
are the same, hence the magnitudes of $1 = $2 and 
$3 = $4 are identical. However, there is no specifica¬ 
tion of the relative phase between superconducting order 
parameters in the two regions. Setting aside more ex¬ 
otic possibilities of phase difference equal to a fraction of 
7T, we are left with two options for the pairing symme¬ 
try: an s-wave, for which the phases of $1 and 4)3 are 
identical, and a d x 2_ y 2 . for which $3 = —$1 (see Fig. 
0 . When only CDW-mediated interaction is considered, 
the two pairing states are degenerate. This has been no¬ 
ticed before^, and it was argued that the degeneracy is 
lifted by other interactions, e.g., antiferromagnetic spin 
fluctuations would favor d-wave. 

We now proceed with the calculation of Tf h . We as¬ 
sume and then verify that the 4>i(fc) = 4>2 (k) = 4>(fc) 
are even functions of momentum k. The linearized gap 
equation (USD. whose solution exists right at T = T C: then 
becomes 

$(fc) =2 TJ2 I Uf(k - k')\G 2 (k')G 5 {-k')^{k'). (20) 

fc' 


In the right hand side of Eq. (12U1) we first integrate over 
momentum transverse to the FS. Neglecting terms small 
in 1 /TV, we obtain 


4>(w m ,fc||) 


g c T \ 7 dk^ 4 ) (w m ,7||) 

V F J 2 7r Wm -* S ( w m. fc ||)l 
1 

X kj + fc[| 2 - 2/3fc||fc|| + 7 c |w m - u ' m \’ 


( 21 ) 


where /3 = (1 — a 2 )/(l + a 2 ). This factor appears in the 
last term in m because fcy and k'^ are parallel compo¬ 
nents of momenta in different segments of the FS, namely 
near hot spots 1 and 2 , respectively. 

A similar gap equation has been analyzed in the con¬ 
text of spin-mediated pairing near SDW QC P n i 46 ' 52 . To 


make this paper self-contained, we present some details 
of the computation of Tf h in our case. 

It is usually more convenient not to solve Eq. ( 1211 ) di¬ 
rectly, but to add to the right hand side of the gap equa¬ 
tion an infinitesimal pairing condensate 4>o and compute 
pairing susceptibility Xpp = 4>/4>o- The transition tem¬ 
perature T c is the one at which pairing susceptibility di¬ 
verges. This approach has an advantage in that the pair¬ 
ing susceptibility can be analyzed within perturbation 
theory. 

The first iteration gives 

*(“"■- T ' 0) = *o(l + A,og^), 

( 22 ) 

where, we remind, ojq c ~ g c is the upper edge of NFL be¬ 
havior. We note that neither the coupling constant g c nor 
1/N directly appear in Eq. (| 22 f) . i.e., once temperature 
is expressed in units of to q c , the renormalization of $0 is 
fully universal. The presence of log 2 term (i.e., one extra 
power of log compared to BCS theory) is the consequence 
of the singular dependence of the fermionic self-energy 
on momentum along the FS. The log 2 term comes from 
momentum range where E(w m ,/c||) ~ iui m / |fc||| [see Eq. 
(HS|]. At g c /v F > |fc||| > ^/ 7 c|w 0 c| the term l/\u' m - 
fc||)| in Eq. (l2ll) scales as |/c||/w m |. To logarith- 
mical accuracy, the momentum integral over fcy in Eq. 
(I32T1) yields | d/cjj/fcjj oc log|w(J, and the frequency 

integral over u' m then yields j T (\og\uj' m \)/\u>' m \duj' m oc 
log 2 (l/T). For the spin-fermion model, this result was 
first obtained in Ref. 11. 

The log 2 T renormalization of the pairing vertex has 
been found in other contexts^ - — . To see how it is 
relevant for T c one has to go beyond one loop order. 
To get the insight, we first consider the “weak cou¬ 
pling” limit by formally replacing the actual coupling 
l = 2a/{a 2 + 1 ) by an effective l e = 2 ea/(a 2 + 1 ) and 
taking the limit e<l. In this limit, series of log 2 renor¬ 
malizations can be summed up explicitly, and the result is 
4> = <I > oe ie ^ 2 ’ r ^ log2 ( U W T ). yyg see that, at the log 2 level, 
the pairing susceptibility does increase with decreasing 
T, but it does not diverges at any finite T. One then has 
to go beyond the log 2 approximation and include sub¬ 
leading O(log) terms. In the weak coupling limit Kl 
this can be done rigorously, along the lines specified in 
Ref. HH, and the result is that subleading O(log) terms 
do give rise to the divergence of the pairing susceptibility 
at a finite Tf h given by 

Tf ~ ujoce- 1/£ . (23) 
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The exponential dependence is the same as in BCS for¬ 
mula, which is not accidental because in the limit e 1 
the main contribution to superconductivity comes from 
fermions away from the hot regions, where self-energy 
has FL form. However, in distinction to BCS formula, 
the prefactor wo c is not the upper cutoff for the attrac¬ 
tion but rather the scale set by the coupling constant g c . 
The proportionality of Xj? h to the coupling g c is the fin¬ 
gerprint of the pairing near a quantum-critical poin t 47 ’ 62 . 

For the physical case e = 1, we expect from C3l) 
T£ h ~ uj 0c . To obtain the exact relation we solved Eq. 
(12T1) directly, using the finite-temperature form of the 
fermionic self-energy. Typical internal momenta and typ¬ 
ical internal frequency are of order, fcjj ~ 7 c u m ~ qwoc 
and ui ~ ujoc . In this situation, fermions from both NFL 
and FL regions contribute to the pairing. For numerical 
evaluation of T£ h we extracted the fermionic dispersion 
in hot regions from ARPES data for B^Si^CaC^Og+rr 
(Ref. HI) and obtained a = 0.074. Using this value for a 


we obtained numerically 

T c ch = 0.0025<? c . (24) 

For comparison, in the spin-fermion model the critical 
temperature at SDW QCP is^ T s c p = 0.0073&, (T c s p ~ 
140 K for g s ~ 1.7 eV). We see that T c in spin-fermion 
and charge-fermion models on top of the corresponding 
QCP are comparable if g c > g s . If this is the case, then 
CDW fluctuations give rise to substantial enhancement 
of the superconducting T c around CDW QCP. One also 
should keep in mind that the result we quoted for T® p 
due to spin fluctuation exchange is T® p right on top of 
SDW QCP. Near the CDW QCP, magnetic £ s is finite 
and spin-mediated T® p is reduced. 

Away from the CDW QCP, charge-fluctuation ex¬ 
change preserves FL behavior and charge-mediated T c 
drops and eventually follows weak coupling BCS formula. 
The self-energy at a finite charge correlation length £ c is 
modified compared to Eq. (USD and is given by 


S(fcj|, cj m ) 


19c [ x 

-sgn(w m ) 

ttvfJc 





(25) 


Now even right at a hot spot (at fen = 0) the self-energy 
has a FL form 

uj 2 

E(w m ,0) = \ c {iu>m) - (26) 

4w c f 

where 


Ac 


gcCc 
2ttvf ’ 


and w c f 



9c a 
4fA 2 a 2 + 1 


(27) 


The dimensionless charge fermion coupling A c (the ratio 
of g c to typical fermionic energy vf Qr 1 ) decreases when 
£ c decreases. Once A c < 1, charge-relaxation scale w c f 
becomes the upper energy cutoff for the pairing, and T c 
follows BCS-Eliashberg-McMillan formula^ 

T c ~ w c t e A <=°. (28) 

We present the numerical result for the behavior of T£ h as 
a function of £ c in Fig. [11] (red line). A similar reduction 
of charge-mediated T£ h is expected on the other side of 
CDW QCP, in the charge-ordered state. 


III. CHARGE-FERMION COUPLING 
CONSTANT FROM THE SPIN-FERMION 
MODEL 


To compare the magnitudes of g c and g s we com¬ 
pute their ratio within a particular microscopic model for 


charge order in the cuprates. Namely, we assume, as in 
earlier works by several groups including ii a 11 ’ 12 ’ 52 ’ 62 ’ 65- — 
that spin fluctations develop at higher energies than 
CDW (and superconducting) fluctuations, and CDW or¬ 
der emerges due to spin-fluctuation exchange, as a com¬ 
posite order. 

To this end, we consider a 2D itinerant electron system 
in which the primary interaction between fermions is me¬ 
diated by soft collective spin fluctuations at the antiferro¬ 
magnetic momentum = (n,n). Such an interaction, 
shown as the wavy line in Fig. [GJ scatters fermions be¬ 
tween hot spots 2 and 4, 1 and 3, etc, and is proportional 
to the dynamical spin susceptibility: 

^eff — U s (q, II m ) ^ ( Cfo a (T OL pCk-\-q,l 3 C!p ^(7'ySCp—q^ 

k,p 

(29) 

where 

Uf{q,n m ) =g 2 s Xs(qMm) 


6 2 + (q- Qtt) 2 + 7 s|Hmf 

where y s is the corresponding Landau damping coeffi¬ 
cient and the scale g s sets the magnitude of T® p for spin- 
mediated superconductivity. 

Except for special cases, there is no rigorously justified 
procedure to obtain Eq. (1^51) starting from a model of 
fermions interacting with some short-range interaction 
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FIG. 6. Schematic representation of spin-mediated interac¬ 
tion. Near an antiferromagnetic quantum-critical point in 
a metal hot fermions scatter into each other by exchang¬ 
ing soft antiferromagnetic spin fluctuations with momentum 
(71", 7r) (the wavy lines). 


U (r) because the main contribution to static part of f/| ff 
comes from fermions with high energies. 

One commonly used approach is to treat static part of 
spin-mediated interaction phenomenologically and just 
postulate Ornstein-Zernike form of the static x s (q 7 0) 
(this is the same procedure that we used in the charge- 
fermion model). Once the model with static spin- 
mediated interaction is established, one can compute the 
dynamical part of Xs(q 7 S2 m ) (the Landau damping term) 
within this model, as it comes from fermions with small 
energies. Within this approach, one cannot relate g s with 
U(r), but one can express Landau damping coefficient 
via g s - The relation is^ = 2g s [2/(7 ruf*)], where for 
comparison with RPA below we pulled out factor of 2 
due to spin summation. 

A complementary approach is to treat t/f ff as the 
charge component of the fully renormalized vertex func¬ 
tion r a7i ^ 5 (q, fl m ) at momentum transfer q near Q^. 
The vertex function T is the opposite of physical anti¬ 
symmetrized interaction (a direct interaction minus the 
one with outgoing fermions interchanged). The vertex 
function can be obtained in RPA by summing up partic¬ 
ular ladder and bubble diagrams which form geometrical 
series (for details see Ref. [69). The approach is best un¬ 
derstood when U (r) is approximated as on-site Hubbard 
interaction U. The RPA gives 


U 


u 


r “^( q ’ - 1 - [/2n 2 (q, n m ) 6ap5lS + 1 - tm(q, n m ) 5a * Sfh 

u „ „ u 


2(1 + £/n(q, fl m )) 




2 (l-£/n(q,fi m )) 


Goi/3 ’ G-y8 


(31) 


where to split the vertex into spin and charge parts we 
used cr a p • = —SapS^s + 2 S a g6p-y For a repulsive inter¬ 

action U > 0, and the interaction in the spin channel is 
enhanced and at large enough U diverges at q, at which 
static n(q, 0) > 0 is at maximum. We assume that the 
maximum of n(q, 0) is at q = Q^. 

Near n(Q,r, 0) = 1/17, the interaction in the spin chan¬ 
nel well exceeds the one in the charge channel, and one 
can keep only the spin component of the interaction, i.e., 
approximate the dressed interaction by Eq. with 

Ufa a p ■ = 2 ( 1 _f/n(q,n m ))^ ' (32) 

Expanding the polarization operator near antiferromag¬ 
netic momentum and zero frequency, we obtain 

n(q,H m ) = n(Q w ,0)-a(q-Q w ) 2 -^^ (33) 

(the last term comes from fermions near the FS and the 
prefactor for term is known exactly). Substituting 


this form into (1321) we obtain after simple manipulations 
the same U° s ( q, fl m ) as in Eq. (ETUI) with 



ZO-7T 


l-un(Q^o) 

uc„ 


Is = 2 g s —r ( 34 ) 




We see that the expression for the Landau damping 
coefficient is exactly the same as in the other (semi- 
phenomenological) approach, the only new element of 
RPA is that g s = 1/(2 C„) is related to the behavior 
of the static polarization bubble. Formally, g s doesn’t 
depend on 17, but in reality C n is of the same order as 
n(Q„0) (in units where lattice constant a = 1), and 
the latter is approximately 1/17 near a SDW QCP. As a 
result, g s is fact is of order U. 

The interaction mediated by spin fluctuations gives 
rise to d-wave superconductivity and to fermionic self¬ 
energy. In the FL regime, E(w m ) ~ A s u m , where 
A s = 3g/(47TOF^J 1 ). This self-energy gives rise to mass 
renormalization m*/m = 1 + A s and to quasiparticle 
residue 1/Z = 1/(1 + A s ). We will include this renor- 
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malization into the calculations below. 

We now proceed to construct the interaction in the 
CDW channel. 

The CDW instability with the ordering momentum 
Qx = {Q, 0 ) and Q y = ( 0 , Q) emerges in this approach as 
a preliminary collective instability at a finite £ s , due to 
spin-fluctuation exchange. A way to obtain CDW insta¬ 
bility is to introduce infinitesimal CDW field Ajj?, which 
couples to incommensurate component of charge density 
as A?4_ -Q/ 2 ,c/a/ 3 c fc+Q/ 2 ,/ 3 ; an d compute susceptibility 
with respect to this field. This has been don o 13 i 14 i 17 
by summing up ladder series of renormalizations due to 
spin-fluctuation exchange. Each act of spin-fluctuation 
exchange transforms hot fermions near, say, hot points 1 
and 2 in Fig. [ 6 ] into another set of hot fermions near the 
points 3 and 4. The set 3,4 is generally different from 
the set 1,2 (because directions of the Fermi velocities are 
different), so to obtain the susceptibility one has to solve 
the set of two coupled equations for fully renormalized 
A^ with the center of mass momentum k either between 
points 1 and 2 or between points 3 and 4. 

There is no rigorous justification why one should re¬ 
strict with only ladder diagrams, even at large N. The 
first non-ladder diagram is small numerically, but not 
parametrically, compared to the ladder diagram of the 
same loop order. Accordingly, there is no point to keep 
N as artificially large parameter, and in this section we 
set the number of pairs of hot spots N to their actual 
value N = 4. 

We present the diagrammatic representation of this set 
of equations in Fig.[7](a,b). The linear “gap” equation for 
A^ has been analyzed in Ref. [lj and the outcome is that 
the CDW susceptibility diverges at a finite T cc j w before 
the system develops CDW order. The critical temper¬ 
ature T cdw decreases as £ s decreases and vanishes at a 
finite critical £ s , setting up a CDW QCP at some dis¬ 
tance away from SDW QCP— (see Fig.[l]). 

Alternatively, one can combine pairs of subsequent 
renormalizations of A 9 into new effective interaction at 
small momentum transfer [see Fig. 0(c)], which we la¬ 
bel as U c and show graphically in Fig. [HJ This com¬ 
posite effective interaction is the convolution of two 
fermionic propagators and two spin-fluctuation propaga¬ 
tors. Explicit calculation show s 21 ' 70 that U c is numeri¬ 
cally smaller but parameter-wise of the same order as a 
single spin-fluctuation propagator - one extra power of 
g s in the numerator gets cancelled out by the Landau 
damping coefficient 7 S oc g s in the denominator. This 
effective interaction is repulsive (U c > 0) because the po¬ 
larization bubble in the particle-hole channel is negative, 
as opposed to the bubble in the particle-particle channel. 

The corresponding term in the Hamiltonian is 

9 U 

2 OctpUfiv T 2°’“/ 3 ' * 

(35) 

where subindices 1 and 2 indicate that the corresponding 




FIG. 7. The linearized “gap” equation for the CDW order pa¬ 
rameter Ajj? ~ ( c !+q/ 2 c *:-Q/ 2 )' We define the center-of-mass 
momentum of hot spots 1 and 2 as fco and that of hot spots 
3 and 4 as k n . Panels (a) and (b): the coupled gap equations 
for and A^ . Panel (c): The gap equation for A^ q only, 
obtained by combining panels (a) and (b). We treat the com¬ 
posite object in the dashed frame as the effective interaction 
U c (the prefactor —1 reflects that interaction appears in the 
diagram with the minus sign). 

momenta are near hot spots 1 and 2 , and the spin factors 
originate from 

(^"7/3 * &as) ( 0^7 * ^’5t') 

7<5 7 )^ ck /3 * ^ 7 < 5 ^ ^ b'-fiii/u ^jji^j 

9 1 

— 2 < W ( W T ' &ilv• (36) 

Only the first, term is relevant to CDW insta¬ 

bility as it renormalizes A^, which acts between fermions 
near hot spots 1 and 2 , with the same spin components 
(i.e., it is convoluted with 5 a p), therefore we can drop 
the 7 a p ■ (Jy. 1 , component in Eq. (l35l) . 

The composite interaction U c can be approximated as 
a constant if the deviation of the fermionic momenta from 
corresponding hot spots (e.g., regions 1 and 2 in Fig. 0 
are smaller than the inverse spin correlation length 
At larger deviations from hot spots, U c becomes function 
of momenta and decreases. 

We now turn to the calculation of the fermionic self¬ 
energy and pairing instability. At a first glance, an in- 


— U c {c\ a C2 t u){c\ ll c\^) 
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1, /3 d y p 4, <5 <r M7 1, /j, 



2, ck 3, 'y ^5i/ 2, zz 


FIG. 8. The diagrammatic representation of the composite 
effective interaction U c (same as in the dashed frame in Fig. 
[3 This effective interaction is a convolution of a particle-hole 
bubble and two antiferromagnetic spin-fluctuation propaga¬ 
tors. 


teraction with a positive (repulsive) U c cannot give rise 
to the pairing instability with sign-preserving gap be¬ 
tween the regions 1 and 2. On a more careful look, 
however, we notice that the interaction in Eq. (1351) is 
the one at small momentum transfer (both incoming and 
outgoing fermions are near the same hot spot), while to 
analyze CDW-mediated pairing (Fig. [5]) and fermionic 
self-energy (Fig.QJc)) one needs density-density interac¬ 
tion at momentum transfer approximately equal to the 
momentum difference between hot spots 1 and 2, namely, 
Q y = (0, Q). To obtain this interaction, we need to inter¬ 
change one creation and one annihilation fermionic oper¬ 
ator. Then we obtain from Eq. (1551) 

He = 2 ~ ( C 2 ,a^O!/ 3 C l,/ 3 )( C l,/i^A“^ C 2 ,y)- (37) 


Viewed this way, the effective interaction is attractive and 
is capable to give rise to pairing. 

Another, more standard way to verify that the density- 
density interaction with momentum transfer Q is at¬ 
tractive is to again extract it from the vertex func¬ 
tion r a#ti/J „(q, Sl m ). To second order in the spin- 
fluctuation propagator, there are two direct and two anti¬ 
symmetrized diagrams for r a(1 ^„ (q, fl m ) which contain 
spin propagators with momenta near Q^. We show them 
in Fig. [9] In evaluating these four diagrams, we addition¬ 
ally require that both spin-fluctuation propagators carry 
the same momenta, as only then one can cancel extra 
power of g s . It is easy to show that only the fourth dia¬ 
gram (the one from the anti-symmetrized part) satisfies 
these conditions. This diagram is exactly the same as the 
one in Fig. [5J but there is an extra minus sign in front 
of it. The evaluation of the diagram itself gives —U c be¬ 
cause particle-hole bubble is negative. The extra (—1) 
in front of this diagram cancels the overall minus sign. 
As a result, the charge component of the vertex function 
becomes 


r c 

x a[A,/3iv 



(38) 


where the spin structure is obtained the same way as in 
Eq. (1361) . Associating the charge component of the vertex 
function with the effective density-density interaction, we 
reproduce Eq. (EZD. 


The effective interaction — (9U c /2)c\ a ci jQ cj M C 2 , M is 
the bare interaction at momentum transfer Q y , and in 
this respect 9U c /2 plays the same role as Hubbard U 
played for our earlier derivation of spin-mediated inter¬ 
action within RPA. Just like we did for spin case, we now 
dress interaction by summing up series of RPA diagrams 
(see Fig. [TUI) . This way we obtain fully renormalized 
(within RPA) effective interaction in the charge channel 


jjeff _ 9E4 _1_ 

c “ 2 l-9E7 c |n c (q,n m )| 


(39) 


Expanding the polarization operator n c near, say Q = 
Q y , we obtain 


|n c ( q ,n m )| = |n(Q„,o)| - c v (q- Q y ) 2 - ln ™ lor + 1 

TTVj? AOi 

(40) 

Substituting this form into (I5U1) we obtain the effective 
charge-mediate interaction the same form as in Eq. ([5]) 
with 


1 ^_ 2 l-9l/ c n(Q„,0) _ = 1 a 2 + l 

9c 2 Cy 9 U c C y ’ 7c 9c ttv 2 f a 

(41) 

To proceed further we approximate the dynamical spin 
susceptibility y s (q, fl m ) by its value at q = Q 7 and 
f l m = 0 and integrate over fermions within the momen¬ 
tum range of the width A around hot spots. Like we said 
above, the approximation of U c by a constant is valid 
when momentum deviations from a hot spot are at most 
of order £~ 1 , so A£ s is generally of order one. We also 
assume for simplicity that CDW order parameter has a 
pure d -wave form, i.e., set our parameter a to be one. 
Within this last approximation, the polarization opera¬ 
tor has the same form between points 1-2 and 3-4. Eval¬ 
uating the polarization operator n c we then find near 
q = QgH 


|n c (q)| 


A 

\Z2it 2 vfQ + A) 


r< (q-Q*) 2 
- y A 2 


(42) 


where C y is of order one and we remind that A = 
3 <? s /(47T?;p’(;“ 1 ) is the mass renormalization due to spin- 
fluctuation exchange. The variable C y , which we intro¬ 
duced in (HUl) is related to C y as 


C v = C v -=— -. 

\[2i r 2 Uj?(l + A)A 


(43) 


Within the same approximation the composite interac¬ 
tion U c is given by 

Uc = ( U 2 ) 2 |n c (Q y )| = f s e k 2 ,, (44) 

V27tWf(1 + A) 

such that 

9[/jn c (Q y) | = (3 + (45) 
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FIG. 9. The vertex function with momentum transfer near Q y in two-loop order (two spin-fluctuation propagators). Other 
two-loop diagrams (not shown) contain spin-fluctuation propagators with small momentum transfer and are irrelevant for our 
purposes. The charge component of this vertex function is 914/2. The opposite of this charge component (i.e., —914/2) is the 
bare interaction in the charge channel at momentum transfer near Q x or Q y . 





FIG. 10. The RPA diagrams for the dressed effective charge 
interaction Each double slid line is the ’’bare” U c - the 
charge component of the vertex function at two-loop order. 
The dressed interaction I/® ff can be viewed as charge fluctua¬ 
tion exchange (see text). 


Using the condition 9f/ c |n c (Q y )| « 1, we eliminate un¬ 
known scale A and obtain 



i/? 


ri _ r< 3g s t; y 

v ~ v 2 ^ 4(1 + 


A ) 2 


(46) 


Hence 

1 _ 3tt 2 /1 + A\ 2 o _ 1 /I + A\ 2 

90 ~ 2Cy ~ 98 16Cy \ A J ^ 9S Cy\ A j 

(47) 

We see that within this approximation the ratio g c /g s 
depends on the value of dimensionless parameter C y . 
To obtain this parameter one needs to know more pre¬ 
cisely system behavior at energies comparable to A. Still, 
if C y (X/(l + A)) 2 < 1, then g c > 2 g s , in which case 
superconducting T£ h from the exchange of near-critical 
charge fluctuations well may exceed T® p from the spin- 
fluctuation exchange. This is the central result of this 
section. 

A more quantitative analysis requires extensive numer¬ 
ical calculations and is beyond the scope of this work. We 
also emphasize that there is no known controllable pro¬ 
cedure of the derivation of the effective interaction medi¬ 
ated by near-critical collective bosonic fluctuations, the 
RPA which we used is an uncontrollable approximation. 
And we also recall that the composite interaction U c (the 
convolution of two fermionic propagators and two spin 
propagators) does depend on external momenta and fre¬ 
quency, and already the calculation of T c( j w requires one 
to solve integral equation for momentum and frequency 
dependent full . 


FIG. 11. The behavior of T/ h ’s as a function of l/£ c with and 
without an external field H, obtained by explicitly solving Eq. 
(BHD - Without a magnetic field T’/ h decreases as becomes fi¬ 
nite and at small enough £ c crosses over from quantum-critical 
to BCS-like behavior [see Eq. (1281) 1. At a finite field, T/ h 
at £c = oo is somewhat reduced, but, most important, T® h 
now vanishes at a finite In numerical calculations we 
used g c = 0.75 Ef, a = 0.076, and [ibH — 5 x 10 ~ a Ef (for 
Ef = 1 eV, this H is ~10 Tesla). For these parameters, 
= 4.35/fcp. The critical is well described by Eq. (1501) 
- plugging £7 into this equation gives 1.2 in the right hand 
side, close enough to the actual 1. 


IV. SHRINKING OF A SUPERCONDUCTING 
DOME IN A MAGNETIC FIELD 


Finally, we discuss the issue of how T)r h , mediated by 
collective degrees of freedom, evolves in the presence of 
an external magnetic field. For definiteness, we focus on 
the role of near-critical charge fluctuations and neglect 
the contribution to T c from spin fluctuations. 

It has been found experimentall y 33 ' 71 that in the pres¬ 
ence of a magnetic field H superconducting T c , viewed as 
a function of doping, splits into two domes, and the one 
at higher doping is centered at or very near CDW QCP. 
As H increases, the maximum of T c in this dome is some¬ 
what reduced, but, most notably, the width of the dome 
shrinks, i.e., superconductivity get progressively confined 
to a CDW QCP. 
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We show that this behavior is reproduced within a 
quantum-critical CDW pairing scenario. The argument 
is rather straightforward right at CDW QCP, T/ h is 
set by charge-fermion coupling g c , and to reduce T£ h 
one would need to apply a rather strong magnetic field 
HbH ~ g c ■ Away from CDW QCP, in the FL regime, 
T£ h is reduced and eventually follows BCS formula. In 
the latter case, a much weaker fi bH is needed to kill 
superconductivity. 

To see how this works in practice, we solved for charge- 
mediated T/ h at a finite charge correlation length £ c by 
assuming that the dominant effect of the field is Zeeman 
splitting of fermionic energies in the particle-particle bub¬ 
ble. Within this approximation, the linearized integral 
equation for the pairing vertex function <I>(w m , &||,) (the 
one which has a solution at T = T^ h ) is 

$(w m ,fe||) 

_g c T r dfcfi $(^m> fc ||) s g n KJ 

v f J 2 tt u'm 

1 

X _ 

k l + k \\ ~ 2 P k \\ k '\\ + 7c|w m -u ' m | + £T 2 
_gcTy^ f dk || <K<o fc ||) l W m -* S ( W m’ fc ||)l 

V F J 2lT [ W m ~ * S ( W m , fc ||)] 2 + ( MB #) 2 

k \\ + k \ 2 ~ 2 P k \\ k \\ +7c|^m - +^ 2 ’ ^ ^ 

where the self-energy is given by Eq. (USD- 

iii the FL regime, when A c = g c /(2nvF^c 1 ) — 1> 
<l>(fc||,w m ) can be, to logarithmical accuracy, approxi¬ 
mated by a constant $, and Eq. (14S1) reduces to 


$ = 


■log 


^cf 


1 + A c 6 (T 2 + H 2 y/ 2 


$. 


(49) 


where u ) c f ~ g c /A 2 has been introduced in E71) . The 
superconducting T£ h becomes zero at a critical A£ r , given 

by 


Ag r 

1 + A- 




log h =1 ’ 


(50) 


or, with logarithmical accuracy, at A“ ~ 
1 /log {g c /g, B H). At smaller A c , i.e., at larger devi¬ 
ations from CDW QCP, there is no charge-mediated 
superconductivity. 

We solved the gap equation numerically and obtained 
T/ h as a function of A c and hbH. We plot the results 
in Fig. 1111 and present the corresponding phase diagram 
schematically in Fig. [l^] We see that, indeed, super¬ 
conducting dome gets sharper in the field, i.e., charge- 
mediated superconductivity gets progressively confined 
to CDW QCP. We did not do calculations on the other 
(ordered) side of CDW QCP, but by generic reason we 
expect a similar shrinking of T c range. The shrinking of 
T c range with increasing held is fully consistent with the 
experimental dat a 33 ' 71 . 



FIG. 12. The variation of the onset temperature of supercon¬ 
ducting order mediated by near-critical charge fluctuations in 
the presence of an external field H. As our numerical results 
show (Fig. 1111) . the range of the superconducting dome shrinks 
as magnetic field increases. 


V. CONCLUSION 


Motivated by the observation of a static charge or¬ 
der in the cuprates and the enhancement of T c at its 
onset, we studied in this work the pairing mediated by 
charge fluctuations around the quantum-critical point to¬ 
wards an incommensurate charge order with momentum 
Qx = (Q,0) or Q y = (0,Q). Our main goal was to un¬ 
derstand whether charge-mediated pairing near a CDW 
QCP yields T c comparable to that obtained from spin- 
fluctuation exchange. 

We first considered a semi-phenomenological charge- 
fermion model in which hot fermions (the ones at the 
FS, connected by Q x or Q y ) interact by exchanging soft 
collective excitations in the charge channel. We obtained 
bosonic and fermionic self-energies in the normal state 
and used them as inputs for the analysis of the quantum- 
critical pairing problem. We found, in agreement with 
earlier works^H that the charge-mediated pairing inter¬ 
action is attractive in both d-wave and s-wave channels. 
The d-wave pairing becomes more favorable once we in¬ 
clude other contributions to the pairing interaction from, 
e.g., antiferromagnetic spin fluctuations. We found that 
the critical temperature T c scales with the charge-fermion 
coupling constant g c , and that fermions from NFL regime 
very near a hot spot and from a FL region further away 
from a hot spot contribute to the pairing. In this respect, 
pairing near a CDW QCP is similar to the pairing by spin 
fluctuations near a SDW QCP. We obtained the value of 
Tc/Sc numerically. 

We next considered the microscopic model in which 
spin fluctuations emerge at higher energies than charge 
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fluctuations and are therefore the primary collective de¬ 
grees of freedom. Charge fluctuations emerge at smaller 
energies as composite fields, made out of pairs of spin 
fluctations. Within this model, we were able to ex¬ 
press charge-fermion coupling g c via the underlying 
spin-fermion coupling g s and relate T£ h due to charge- 
fluctuations near a CDW QCP to T® p due to spin fluctu¬ 
ations. We found that, at least within RPA, T£ h due to 
charge fluctuations is comparable to T® p that due to spin 
fluctations and may even exceed it, i.e., superconduct¬ 
ing T c does get a substantial enhancement near a CDW 
QCP. 

Finally, we analyzed the behavior of charge-mediated 
T c in the presence of a magnetic held and found that 
the dome of Tf 1 around a CDW QCP indeed shrinks 
as magnetic held increases, because a held destroys 
superconductivity faster in non-critical regime than in 


the quantum-critical regime and hence enhances charge- 
fluctuation component of T c near a CDW QCP. 

This result and the one that the contribution to T c 
from critical charge fluctuations can be larger than the 
contribution to T c from non-critical spin-fluctuations, de¬ 
spite that charge fluctuations are by themselves made 
out of spin fluctuations, may explain the experimental 
observation that in a magnetic held T c gets progressively 
confined to the doping range around the doping at which 
charge order likely emerges at T = 0. 
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